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Abstract. 

The cyclic evolutions and associated geometric phases induced by time-independent 
Hamiltonians are studied for the case when the evolution operator becomes the identity 
(those processes are called evolution loops). We make a detailed treatment of systems 
having equally-spaced energy levels. Special emphasis is made on the potentials which 
have the same spectrum as the harmonic oscillator potential (the generalized oscillator 
potentials) and on their recently found coherent states. 
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Since the appearance of Berry's work [1], much effort has been spent in studying 
geometric aspects of nonrelativistic quantum mechanics [2 — 8] . In particular, to any cychc 
evolution of the vector state, \iI^{t)) — e^^\ip{0)), there has been associated a geometric 
phase 



where r is the period of \i^{t)) e H, (^/'(t) = 1, (/» e R, is the Hilbert space of 

vector states of the system, and H{t) is the Hamiltonian [3]. /3 describes global curvature 
effects arising on the space of physical states, which is the projective space V formed by 
the rays or the density operators instead of Ti. Due to this curvature, the horizontal 

lifting (parallel transport) of the closed trajectory \i/j{t)){'i/j{t)\ G V leads to a trajectory 
\ipH{t)) which is, in general, open on 7i. The holonomy of this lifting is the Aharonov- 
Anandan geometric phase factor e*^. 

The phase determined up to a multiple of 27r, generalizes the Berry phase, which 
originally was defined just for adiabatic cyclic evolutions [1] . Further generalizations of P 
have been designed and can be found in the literature [3,5 — 6]. There is also a lot of 
work dealing with the calculation of the geometric phases when the Hamiltonian is time 
dependent (either explicitly or implicitly through certain sets of time-dependent param- 
eters [9 — 17]). In this paper, we will address the study of the geometric phases when 
the Hamiltonian is time independent, i.e. H{t) = H. This choice is done because, it 
seems to us, there is a widespread belief that the geometric phases appear only when the 
Hamiltonian is time dependent, which is wrong. This is, perhaps, motivated by the his- 
torical development of the subject and the following reasoning: the eigenstates \En) of H 
evolve according to \En{t)) = e~^^'^*/^\En), where the En are the energy eigenvalues and 
n denotes a set of discrete subscripts. These evolutions are cyclic with period (arbitrary) 
T and (j) = —EnT/h. Therefore, from (1), P = for these states, and as usually the only 
cyclic states at hand for these systems are the eigenstates of the Hamiltonian, one is led 
to the wrong conclusion stated above. However, Aharonov and Anandan found nonnuU 
geometric phases for a spin 1/2 in a constant homogeneous magnetic field [3]. The same 
will be true for any other two-level system described by a time-independent Hamiltonian 
[11, 15]. For nonspin systems it is possible to prove the existence of nontrivial geometric 
phases for the harmonic oscillator [15, 18] and some other physically interesting models 
(such as the localized states of an electron on a crystal [15]). 




(1) 
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For independent reasons, in order to be used as the starting point for the techniques 
of "controhng" and "manipulating" the quantum systems, the evolution loops (EL) were 
proposed (although without this name at that time) in 1977 by Mielnik [19] and further 
developed by him in 1986 [20]. Those loops are specific dynamical processes induced either 
by time-dependent [19 — 21] or time-independent [22] Hamiltonians, for which the evolution 
operator U{t) becomes the identity 1 (modulo a phase) for a certain time r > 0, i.e.: 

U{t) = e'H, (2) 

where U{0) = 1 (see also [23]). The EL are useful because when perturbed by some 
additional external fields, the system can be driven to attain any desired unitary operation 
on H due to the accumulation of the small precessions of the distorted loop [19 — 21]. In 
the context of geometric phases a system performing an evolution loop is very interesting 
because any state becomes cyclic at t = r: 

|V'(r)) = e^'^|V'(0)). (3) 

Therefore, it could (we will show) have an associated nonnuU geometric phase. We will 
restrict ourselves in this paper to the evaluation of the geometric phases associated to an 
evolution loop when the Hamiltonian is time independent. 

Suppose one has a system with a time-independent Hamiltonian H whose evolution 
operator performs an evolution loop. Hence, any vector state \ip{t)) comes back to itself 
at t = T (see equations (2-3)), and its geometric phase can be easily evaluated because the 
evolution operator U (t) = e"*^*/'^ commutes with H: 

/5 = + n-'/ {ii;{0)\U^{t)HU{t)\i;{0))dt = (l) + n-\{H), (4) 
Jo 

where (H) = (V'(O)liy |V'(0)). Expressing |V'(0)) in terms of the basis {\Em)}, \ip{0)) = 
Y^^Cm\Em), with Cm = (-E^|V'(0)), wc find that (4) becomes: 

P=(j) + n-\J2\'^rnf^m. (5) 

m 

Note that formulas (4-5) are applicable to the cyclic evolution of a vector state induced 
by any time-independent Hamiltonian regardless of whether or not the system performs 
an evolution loop. However, if the system has an evolution loop, then (4-5) will be valid 
for any initial condition. In particular, for |V'(0)) = \En), i.e., Cm = Snm, it turns out that 
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(j) = —EnT/h, and hence = 0. If at least two c^'s are distinct from zero, however, the P 
associated to the corresponding cychc state will be, in general, nontrivial (see [11], section 
3.1). 

There are in the literature some interesting systems whose time-independent Hamil- 
tonian induces evolution loops [15, 18,22 — 23]. Here, we will show the existence of an 
evolution loop for Hamiltonians whose spectrum consists of equally spaced energy levels 
of the form: 

En = Eo+ nAE, (6) 

where AE > is the constant spacing between the levels and Eq is the ground state energy. 
The subscript n e Z+ takes values in [0, AT], where N is finite if H is finite-dimensional or 
infinite if H is infinite-dimensional. The evolution operator of this system reads: 

N 

Uit) = Y,e-'^-'/'^\Er.){Er,\. (7) 

n=0 

It is easy to see that the evolution loop is present at r = 2'nfi/ AE: 

N 
n=0 

By comparing with (2), we obtain = —2'kEq/AE. Moreover, according to (4-5), the 
geometric phase for the cyclic state \i^{t)) is: 

/? = 27rii^^^ = 27rf;n|c.P>0. (9) 

n=l 

Notice that the component cq of |^/'(0)) along the ground state \Eq) is not explicitly present 
in (9). If 13 is restricted (modulo 27r) to the interval [0, 27r), then equation (9) admits the 
following interpretation: /3 measures the "energy excess" (in dimensionless units) of {H) 
above the nearest lower energy level E]^ (see Figure 1). If E]^ is given and is changed 

so that Ek < {H) < Ek+i, then to the end /3 = corresponds cyclic states with (H) = Ek 
(this includes in particular \ip{0)) = \Ek))- To any other (3 G (0, 27r) corresponds cyclic 
states with (H) / E^ and vice versa (here necessarily \ip{0)) ^ \Ek))- 

One of the interesting systems with equally spaced energy eigenvalues for which our 
treatment can be applied is a spin j interacting with a constant homogeneous magnetic 
field B, where j > can be either integer or half-integer. Suppose, for simplicity, that 
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the magnetic field points in the z direction, B = Sk. The spin Hamiltonian can be 
expressed as H = —/iJ ■ B = — cUcJs, where fx is the spin magnetic moment, J is the 
spin operator whose components satisfy [Jfe, J/] = ihekin^n, with k,l,n = 1,2,3, and 
liB = a;c > is the precession frequency of the spin around k. We work in the basis 
m) : JsIj, = m/i|j, m), —j<m< j}. Hence, dim('H) — 2j + 1 = N + 1. Due 
to the minus sign in the Hamiltonian, the identifications AE = fiLJc, Eq = —jhuc, and 
\En) = \ j,j — n) with < n < 2j = N are consistent with equations (6-9). Therefore, the 
system performs an evolution loop and so any spin state evolves in a cyclic way, with an 
associated geometric phase given by equation (9). This is true, in particular, for the spin 
j = 1/2. In this case, it has become a convention to express the generic initial state in 
terms of the spherical angles 9,(fi: 

1^(0)) = e-"^/^ cos(^/2)|l/2, 1/2) + e"^/^ sin(^/2)|l/2, -1/2). 



As the ground state in our notation is \Eo) = 
to the geometric phase is ci = e**^/^ sin(^/2). 
phase becomes the usual one [3]: 



1/2, 1/2), the only coefficient contributing 
Therefore, by applying (9), the geometric 



P = 27r|cip = 27rsin2(^/2) = n{l - cos d). 

Moreover, as is well known, a problem involving just two energy levels can be treated as 
a fictitious spin 1/2 interacting with a homogeneous magnetic field [11, 15]; hence, taking 
care in making a judicious identification of the parameters, the same formulas to evaluate 
its geometric phases may be applied. 

At this point, it is worth discussing a geometric interpretation applicable to systems 
with energy levels given by (6). It is easily understood for the spin-1/2 system of the 
previous example, for which the space of physical states (the projective space) coincides 
with the unit sphere S'^ on R^. Any spin-1/2 state is precessing around the z axis, 
performing cyclic evolutions with a geometric phase which is, in general, distinct from 
zero. There are two states, however, for which the evolution is trivial: during the course 
of time they remain static at the north and south poles on S'^, and correspond to the spin 
aligned along and in the opposite direction of the magnetic field. The geometric phase for 
both of them is zero. For a system with N + 1 equally-spaced energy levels, however, we 
have at hand a more interesting (and more complicated) situation: now, instead of having 
two static points on S'^ there are N + 1 points remaining static under the evolution on V 
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(those associated to the with N either finite of infinite). For each one of them the 

geometric phase is zero. Any other state will move across those points performing a more 
complicated (cyclic) evolution on V with a nonnuU geometric phase (in general) which can 
be easily evaluated using equation (9). 

We proceed now to the analysis of another system having equally-spaced energy spec- 
trum. It can be called the generalized oscillator (GO) because its energy levels are exactly 
the same as the ones of the harmonic oscillator Hamiltonian. The GO potentials were 
discovered by Abraham and Moses [24] using the Gelfand-Levitan formalism [25] and were 
generated by Mielnik [26] using a generalization of the well-known factorization method 
[27] (see also [28]). Because of its didactic value, we will point out some steps used in the 
generalized factorization to generate the GO potentials. We will work from now on in the 
coordinate representation with dimensionless units fi = m = u = 1. 

The classical factorization method applied to the oscillator consists in expressing the 
Hamiltonian 

as the two products 

aa^^H^]-, a^a^H-^, (11) 

where a and are the ordinary ladder operators a = {l/y/2){d/dx + x), = 
{l/V2){—d/dx + x) with [a, a^] = 1. The eigenfunctions and eigenvalues of the harmonic 
oscillator can be constructed using the relations 

Ha^ = a^(H + 1), Ha = a(H-l). (12) 

There is a normalized ground state ipoi^) with eigenvalue Eq = 1/2 which satisfies a'ipo{x) = 
=^ i(^o{x) oc while the normalized eigenfunction ipn{x) associated to the eigenvalue 

En — n + 1/2 is related to the ground state through: 

= ^V'o(x). (13) 
vnl 



The generalized factorization method [26] consists in looking for more general operators 



satisfying just the first one of relations (11): 

bb^ = H + ^. (15) 
Hence, the unknown function P{x) obeys the Riccati equation 

/3' + = 1 + a;^ (16) 

whose general solution is 

Pix) = x+ , AeR. (17) 

A + Jo e y dy 

Now, the point is that the product b'^b is no longer related to the harmonic oscillator 
Hamiltonian, but it leads to a new operator Hx: 

b% = Hx-\, (18) 

where 
with 

The requirement |A| > assures that V\{x) has no singularities. The relationships 

analogous to (12) provide the way to obtain the eigenfunctions and eigenvalues of Hy. 

Hxb^ = b\H+l), Hb = b{Hx-l). (21) 

Hence, the states 6n{x) = b^il)n-i{x)/ ^/n^ n = 1, 2, • • •, are orthonormalized eigenfunctions 
of H\ with eigenvalues En = n + 1/2. However, the set {9n{x)^ n= 1, 2, • • •} is not yet a 
basis of L^(R). There is a missing unit vector 9o{x) which is orthogonal to all the vectors 
^^(a;), n = 1, 2, • • •. It turns out to be an eigenfunction of H\ with eigenvalue Eq = 1/2 
satisfying b9o{x) = 0, and taking the form: 

^o(a;)ocexp(^-^ /3(y)ciy^ • (22) 

As the set {9n{x), n = 0, 1,2, • • •} forms a basis in L^(R), then {Hx : |A| > \^/2} is a 
family of Hamiltonians distinct from the harmonic oscillator one but which has exactly 
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the same spectrum as the oscillator has. In the limit |A| — > oo, the harmonic oscillator 
potential is recovered, Vx{x) when |A| oo. 

All the relationships involving the evolution loops and the geometric phase (equations 
(6-9)) can be applied to the GO Hamiltonian (19-20) with Eq = 1/2, AE = 1, r = 27r, (f) = 
—TV, and = oo. In particular, the geometric phase is = 2n{{Hx) — 1/2), and when 
applied to the wavefunctions of the basis {9n{x),n — 0, 1, 2 ■ ■ ■} we recover again (3 — 2n7r. 
Is there any other set of generic states of the GO potential for which we can evaluate 
explicitly the geometric phase? The answer turns out to be positive when considering the 
family of recently found coherent states for the GO Hamiltonian [29] . Here, we will present 
some details of its derivation (for work involving coherent states and their geometric phases 
see [9,11,15,18,30- 32]). 

In the construction of the coherent states of Hx, denoted as \z) with 2; e C, we need 
to identify the "annihilation" and "creation" operators of the system. Because bOn{x) oc 
il^n-iix) ^ ab9n{x) oc •0jj_2(a;) =^ b^ab9n{x) oc 9n-i{x), and an obvious choice is: 

A = b^ab, = b^a^b. (23) 

The coherent states can be defined now as the eigenstates of the annihilation operator 
A with eigenvalues z, i.e. A\z) = z\z). Expressing 1^) in terms of the basis {|^n), n — 
0, 1,2---}, and substituting explicitly that expression in the previous one, we find the 
following family of coherent states (after normalization): 

Vo^2(l, 2; 1^1^) ^ n\^y{n + 1)! 

where \dn) is the ket representing the eigenfunction dn{x) and 0-^2(1, 2;j/) is a generalized 
hypergeometric function defined by: 

^-^ r{a + n)r(p -I- n) n\ 



n=0 



with r(-) the Gamma function. To each value z corresponds one and only one co- 
herent state. However, 2; = is a doubly degenerate eigenvalue of A with two orthogonal 
eigenvectors which will be denoted |^o) and I2; = 0) = By choosing an appropiate 

measure in the complex plane, it can be shown that the set {|^o)) l-^)} is complete in H. 



The relationships presented so far are sufficient for our purpose of evaluating the 
geometric phase. To this end, we need to find the expected value of Hx in the state \z). 
A direct calculation leads to: 

(Hx) = {z\Hx\z) = 1/2 + (26) 

0^2[i-, ^, \Z\ ) 

Finally, substituting (26) in the equation for we obtain the following expression for the 
geometric phase Pgcs of the generalized coherent state: 



Pgcs = 27r 



oi^2(i,i; 


\z\ 




oi^2(l,2; 


\z\ 


P) 



(27) 



To have an idea of the behaviour of Pgcs, we plot it versus Re{z) x lni{z) in Figure 
2. As we can see, the geometric phase is independent of A and depends on 2; in a quite 
different way compared with that of a standard coherent state (SCS) of the harmonic 
oscillator, for which Pscs = '2t:\z\'^ [15,18] (see also Figure 2). This occurs because the 
generalized coherent states discussed in [29] do not tend to the standard ones when A — > 00 
even though the generalized potential tends to the harmonics oscillator potential in this 
limit. A deeper analysis shows that the difference rests on the fact that in this limit the 
annihilation operator A^^ = lim^^oo A = a^a^ is distinct from the standard one a. The 
generalized coherent states, however, could be useful in future applications because the 
product of the uncertainty of the X and P operators for these states is almost minimum in 
this limit 1/2 < lim|;s^|^oo AXAP < 3/2. The question of whether or not there is a family 
of coherent states of Hx tending to the standard ones when A ^ 00, the geometric phases 
included, is open. 

The author acknowledges CONACYT (Mexico) for financial support. 
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Fig.l Schematic representation of the N + 1 energy levels for a system with equally- 
spaced spectrum. If the geometric phase /? is restricted to the interval [0, 27r), then it can 
be interpreted as an energy excess of the system with respect to its nearest energy level 
Ek (the nearest below) in dimensionless units. 

Fig. 2 The geometric phases associated to the standard coherent states of the harmonic 
oscillator (Pscs) and the coherent states of the generalized oscillator (Pgcs) as functions of 
the complex variable z. The minimum values of Pscs and Pgcs are and 27r respectively, 
both at 2; = 0. The missing sections in both surfaces were removed to show the behaviour 
close to the minimum. 
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